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NOMENCLATURE 


Most of the symbols listed here are more fully defined in the text. 
Unless otherwise indicated in the text, the symbols have the units 
listed here: 


Symbols 


Units 

C D 

drag coefficient 


c 

e 

ballistic entry parameter 

ft 3 /lb 

D 

diameter of canopy, length of shroud line 

ft 

f 

breaking tenacity of fabric 

gm / denier 

g 

gravity acceleration 

ft/sec 2 

H 

atmospheric density scale height 

ft 

J 

mechanical equivalent of heat 

ft-lb/BTU 

N 

s 

number of shroud lines 


p 

max 

maximum dynamic pressure 

lb /ft 2 

R 

radius 

ft 

Q 

heat flux 

BTU/ft 2 se 

T 

absolute temperature 

deg R 

V 

velocity 

ft /sec 

lv| 

deceleration 

ft/sec 2 

W 

weight 

lb 

x 

coordinate along symmetry axis 

ft 

y 

coordinate normal to symmetry axis 

ft 

z 

altitude 

ft 

z 

max 

altitude of peak dynamic pressure 

ft 
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radian 


a 


P 

Y 

6 


X 

X 

c 

M- 

P 

P c 

Pt 


M 


Subscripts 


a 

bl 


c 

e 


angle between incident Newtonian 
flow and surface normal 

ballistic coefficient 

angle between velocity vector and local 
horizontal 

effective thickness of canopy 

gray-body emissivity 

molecular mean free path 

molecular mean free path in air at 
sea level 

fractional weight 

density 

density of air at sea level 
density parameter of Earth atmosphere 
density parameter of Mars atmosphere 
str es s 

Stefan-Boltzmann constant 
axial 

boundary layer 

canopy 

pntrv 

' - J 


XT' 


j \ 

r,dr in 


fm free molecule flow 

h hoop 


slugs /ft 2 
radian 


ft 


ft 

ft 


lb /ft 3 
lb /ft 3 
lb /ft 3 
lb /ft 3 
lb /ft 2 

BTU/ft 2 sec deg 4 
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M 

Mars 

max 

maximum 

min 

minimum 

r 

ring 

s 

shroud 

tr 

transitional flow 


The abbreviation log refers to natural logarithm. 
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1. SUMMARY 


The feasibility of a parachute entry system for a Mars lander has 
been investigated in accordance with a work statement appended to 
this report. The principal restrictions in the work statement specify 
the range of ballistic coefficients 10 > (3 slugs/'ft > 10 , the 

weight of the deceleration system limited to not more than payload 
weight W_ = 50 lb , entry velocities in the range 6 < kft/sec <16. 
Furthermore, consideration of a high-temperature high-strength 
plastic (Nomex) undergoing "laminar" peak heating is specified in 
the work statement. 

It was found that with the above restrictions, feasibility limit condi- 
tions are reached for entry velocity V = 6. 3 kft/sec in a parachute 

e -2 .2 

system with ballistic coefficient (3 cs 2x10 slugs /ft . At this feasi- 
bility limit most of the system weight is needed to provide a "massive" 
blunted ieading edge for canopy, so as to alleviate the laminar heat 
transfer, and thereby prevent the development of temperature in excess 
of the permissible maximum for Nomex, T ^ 500°F, at the leading 
edge. 

It was found, further, that greater entry speeds are feasible if the 
restriction to "laminar" peak heat transfer is relaxed and design 
configurations with sharp leading or windward elements are taken 
into consideration. Therefore, additional configurations with peak 
heating in the free molecule flow and/or transitional flow regimes 
have been also analyzed in a first approximation. In these configura- 
tions the principal design parameters, canopy leading edge radius R^ 
and shroud line cross sectional radius R are constrained bv two 
conflicting requirements: (i) they must be sufficiently small to assure 

the occurrence of peak heat flux in the free molecule flow regime and 
(2) they must be sufficiently large to provide structural strength against 
the peak deceleration loads. 
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Both requirements can be satisfied with all-Nomex parachutes over 

the range of (3 considered here. Furthermore, at the lowest value 

P = 10 3 , the feasible maximum entry velocity is V e = 9.7 kft/sec, 

a significant increase over that found for the laminar heating pulse. 

-3 

The (3 = 10 configuration is feasible, for example, with the follow- 
ing design dimensions; canopy thickness 5 and leading edge radius R r 
equal to about . 001 in; shroud line radius R g = . 0036 in, number of 
shroud lines = 3 00; diameter of canopy, taken equal to the length 
of the shroud lines, is D = 76.5 ft. The corresponding canopy weight 
W c = 23.3 lb and shroud weight = 0. 6 lb total less than half of 

the maximum chute weight of 50 lb under consideration. It is to be 

-3 

noted that, for |3 = 1 0 , the theoretically feasible minimum dimen- 

sions 6, R g , etc. , adequate for peak dynamic load, require a minimum 
Nomex parachute weight of about 1 lb. 

The low-weight requirements of the design examples summarized 
suggest that, by utilizing nigh-temperature, high-strength metais 
for parachutes, having the sharp windward edge dimensions to 
assure peak heat load in free molecule flow, the feasible entry 
velocities may be significantly increased without exceeding the 
system weight limitation of 50 lb. The possible extension of the 
velocity limit without exceeding the weight limit is illustrated in 
the following table which compares Nomex (T max = 500°F) perform- 
ance with that of steel (T = 1000°F). 

max 

Maximum with Heating Pulse in Free Molecule Flow 


(3 slugs/ft 2 

io _1 

10' 2 

io ~ 3 

V ft/ sec (Nomex at 500°F) 
e 

2090 

4510 

9700 

V g ft/sec (Steel at 1 000°F) 

3660 

7890 

17,000 


2 


NATIONAL ENGINEERING SCIENCE COMPANY 
S-424 


2. INTRODUCTION 


This report presents the results of an analytic investigation performed 
by National Engineering Science Company (NESCO) under Jet Propulsion 
Laboratory Purchase Order No. EC437910. The contract task require- 
ments are specified in a work statement included here in the Appendix. 

As indicated in the work statement, the objective of the investigation is 
to perform a first order parametric analysis which defines the feasibility 
limits of a Ivlars entry capsule decelerated by an annular parachute. 


The assumptions employed are summarized in Section 3. These assump- 
tions serve to define, in the first place, the general configuration of the 
entry system, and the structural elements which are subjected to critical 
peak dynamic and heat loads during entry. Furthermore, the assumptions 
specify the simplified deceleration and heat transfer laws employed in the 
analysis, consistent with the scope of the first order analysis. Justifica- 


tioi 


iptions is detailed in the text 'where the development 


of the simplified entry equations is presented. 


Starting with the second order shallow entry equations, a description of 
the trajectory is obtained in Section 4 expressing the relation between 
speed V and altitude z in an elementary form convenient for the 
analysis of peak dynamic loads and heat loads. The critical stress for- 
mulas associated with peak deceleration loads are given in Section 5. 

The peak heat loads are determined in Section 6 for the various regimes 
of free molecule flow, transitional and continuum flow encountered in 
descent through the Mars atmosphere. The feasibility limits, illustrated 
for Nomex in particular, are determined in Section 7, by expressing the 
peak temperature as function of entry velocity for a range of ballistic 
coefficients and deceleration system weights under consideration. The 
results are discussed in Section 8, followed by remarks on shock impinge- 
ments in Section 9. 

The salient conclusions, responsive to the work statement in the Appendix, 
are summarized in Section 10. 
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3. BASIC ASSUMPTIONS 


3 . 1 Configuration and Major System Components 
The usual definition of the ballistic coefficient 


p = 


w/g 


(3.1) 


may be converted into the following expression for calculating the range 
of diameters D of interest in this investigation: 


D = 


1 it i ir (3 


(3.2) 


where C_ = 1/2 has been substituted. Since W may vary over the 
range 

50 lbs = W < W < 2W = 100 lbs (3.3) 

p ~ p 

-3 - 1 2 

and (3 ranges from 10 to 1 0 slugs/ft , we obtain the following 
extreme values of D : 


D . 

min 


D 

max 


VI 


8 (50/32. 2) 
10 _1 


A T, ft 


4 


8 gSIl = 89 
11 10" 3 


(3.4) 


In view of the relatively large diameters D < D available for low-0 

~ max K 

entry systems, the drag on the conical canopy is generally incurred 
under continuum flow conditions even at extreme altitudes in the low- 
density Martian atmosphere. Accordingly, by considering that the 
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coefficient =1/2 is associated with Newtonian continuum drag, 

it is inferred that the angle a between the incident axisymmetric 
wind and the canopy surface is a = 45° . 

Three major components are envisaged in the configuration of the 
deceleration system: (a) shroud lines of length D , attached to the 

(b) conical canopy extending between the leading diameter D and 
trailing diameter D 1 = D/\/2 , in view of the given projected canopy 

1 TT 2 

area of D ) , and (c) a structural ring or torus of diameter D , 

with cross sectional radius R , for attachment of the shroud lines to 

r 

the canopy. A sketch of the system is shown in Fig. 1. 

As we shall see from heat transfer considerations, it may be desirable 

to reduce the cross sectional radius R to a minimum value, in 

r 

order to realize free molecule flow heat transfer at the higher alti- 
tudes where peak heat loads develop. Thus, if direct attachment of 
the shroud lines to the canopy can be achieved, the torus may be 
altogether absent, leaving the leading edge of the canopy relatively 
sharp, with cross sectional radius comparable with the canopy thick- 
ness. The feasibility limit of entry at higher velocities V is 
crucially dependent on the possibility of realizing low-(3 systems with 
direct attachment of fine shroud lines to a thin canopy. 

3 . 2 Trajectory Equations 

The entry trajectory is described by the so-called second order solution 
(Ref. 2) in the approximation of small variations of the velocity vector 
relative to the local horizontal. In this approximation we assume in 
effect that the velocity is dependent primarily on the local altitude in 
accordance with the Gaziey formula for constant ballistic coefficient 
bodies descending in an exponential atmosphere. 
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CANOPY 
RING (TORUS) 


SHROUD LINE 


Figure 1 

Sketch of annular parachute, showing general configuration 
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3 . 3 Dynamic Loads 


It is assumed that the critical dynamic loads arise at the instant of 
maximum deceleration. The critical stresses of interest are those 
arising in (1) the leading edge of the canopy and (2) the shroud lines 
at their attachment to the canopy. 

3 . 4 Heat Transfer 

In regard to heat transfer in the free molecule flow regime, it is 
assumed that all surfaces have an accommodation factor of unity and 
that molecular impact is followed by diffuse reflection. 

In regard to continuum heat transfer, it is assumed that stagnation 
heat transfer formulas based on air data can be applied to the Mars 
atmosphere. 

A special interpolation formula between free molecule and stagnation 
heat transfer is postulated for transitional flow heat transfer. 

3 . 5 Maximum Temperature 

It is assumed that maximum temperatures are attained under quasi- 
equilibrium between aerodynamic heat influx and gray-body thermal 
re-rauiation. Uniform local temperatures are assumed to prevail 
throughout any cross section of the shroud lines or of the canopy lead- 
ing edge (ring or torus). Gray-body emissivities close to unity are 
assumed in numerical applications. 
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4. SIMPLIFIED TRAJECTORY EQUATIONS 


The second order solution (Ref. 2) pertinent to small angles of 
inclination y is employed in the calculation of entry trajectories. 
For exponential density variation of the type 


- z/h m 

P =p M 6 


(4.1) 


the velocity V and inclination y are interrelated by 


V = V exp - c 
e e 


(y - y ) sin y 

(p - p ) 

cos y - cos y I e 


1 


cos y 


cos Y h /g R 

M / & M M 


/ g M R M ,\A p e\ 

r mV~“P A 1 p/ 


(4.2) 


(4.3) 


where 


c 

e 



2(3 sin y e 


(4.4) 


The parameters of the exponential atmosphere introduced in Eq. (4.1) 
are determined from density data shown in Fig. 2 for the VM-3 and 
VM-7 engineering models. The VM-3 model data shown in Fig. 2 for the 
z range exceeding 200 kft is based on trajectory computer data furnished 
by JPL. The VM-7 density is about one-half that of the VM-3 in the 
altitude range shown. 


* JPL Trajectory Data Run 109 for R. Weaver, 8/27/66, communicated 
by E. Laumann. 
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Table of Mars Atmosphere Parameters 


Parameter 

H kft 
M 

1 ° 5 P M Slug 

VM-3 

50. 5 

4. 90 

VM-7 

50. 5 

2. 45 


It is interesting to note that, notwithstanding the low surface density 

-3 3 

parameter p^. compared with that on earth, p^_, = 3.4x10 slugs /ft 
(Ref. 1, p. 13), the high altitude Mars atmosphere is favorable to 
deceleration effects owing to the larger Mars scale > H^, = 22 kft. 

The data on the above table, in fact, imply that the Mars density exceeds 
that of the earth at altitudes above 165 kft for the VM-3 and 192 kft for 
the VM-7 models, respectively. 


Considerable simplifications are effected in Eqs. (4-2, 4-3) by noting 
from Eq. (4.3) that under representative entry conditions of interest, 
the variation of cos y in the high altitude trajectory is small down to 
the levels of peak heating and peak deceleration. This is readily seen 
by considering Mars data 

gj^.= 12.3 ft / sec ^ R^ = 2120 mi V^M^M ^ esc * ve ^* ) = 1 .66 • 1 0^ ft /sec 


3 J 1 4 3 ^3 „ 

ctnu. utJteicicucu apccua 


Iv 


-J00O f 4 ■ / c 

v v w i.b/ o oo 


.veil below 


+v>< 


critics.! 


stress and heat loading conditions of interest. For the preceding data 
Eq. (4.3) yields 


The lowest entry speed of interest V = 6000 ft/sec is employed 
here, since for larger V , the inclination y differs from y b 
less than the extreme difference shown in Eq, (4. 5). ' 


10 


I lUfNML. L.l’«UII’«L.un!INU WtVirjHify I 


S-424 


cos y = 


(4.5) 


cos y 


1 + 


9.55 /(1.17) 2 A 

2120 l 3 2 I 


~ 0. 94 cos y 


We therefore find that for the shallow entry angle = 15 the bracketed 
quantity in Eq. (4. 2) is very close to unity in the altitude and speed 
regime of interest. Thus, we find from Eq. (4.5) 


(Y ' V J sin Y p 

92 < -2 — < 1 

— cos y " cos Y — 


(4.6) 


In view of the preceding relation, the velocity equation (4. 2) reduces to 
the familiar expression obtained by Gazley (Ref. 3) for steep reentry: 


V = V exp 
e 


C e P 


(4.2') 


where we have set p << p in the deceleration regime of interest. We 
shall employ Eq. (4. 2’) in the calculation of peak stresses and heat loads 
which, as has been shown, depend primarily on the altitude z implicit 
in p and only weakly on the local inclination y • f° r accuracy in 
trajectory calculations, beyond the scope of the present study, both z 
and y would have to be determined by simultaneous solution of 
Eqs. (4. 2) and (4. 3). 

A special consequence of the simplified trajectory in Eq. (4. 2') is that 
the maximum deceleration, associated with maximum dynamic pressure 
pV 2 /2 , is independent of P fCf. Eq. (3-20) of Ref. l] : 


1 '.1 

dV 

V 

| |max 

dt 


V sin y 
e ' e 


max 


2 e H 


M 


(4.7) 
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However, the altitude z = z at which the peak deceleration occurs 

max 

depends on (3 as given by [Cf. Eq. (3-30) of Ref. l] 


z 

max 



log 


H M P M 
P sin y e 


(4.8) 


4 

For example, for representative data of interest V = 10 ft/sec, 
sin y = 0. 259, (3 = 0. 01 slugs/ft" 1 , we find the peak deceleration of 

92. 0 ft/sec occurring at the altitude 346 kft in the VM-3 and an 
altitude of 3 1 2 kft in the VM-7 atmosphere models. 


We conclude this section by comparing the V versus z trajectory based 
on Eq. (4. 2') with the JPL 6-degree of freedom trajectory result for 
(3 = 10 ^ at entry conditions Y e = 20° , = 13, 190 ft/sec , 

z^ = 780 kft. The agreement in the comparison shown in Fig. 3 
is plausible, indicating that for the purpose of the present calculations 
the analytic results of this section will be sufficiently accurate. The 
analytic results provide, of course, the advantages of utmost simplicity 
in the calculation of maximum heat load and stress load during entry. 
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5000 


V ft/ sec 


10,000 


15,000 


Figure 3 

Comparison of calculated entry trajectories (£ = 10" 2 slugs/ft 2 , 
7 e = 20 deg, z e = 780 kft, V = 13. 19 kft/sec, 

VM-3 atmosphere model) 
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5. MAXIMUM DECELERATION STRESSES 


In this section we describe the maximum stresses associated with 
maximum deceleration loads on the entry system sketched in Fig. 1. 

In particular, we consider the hoop stress in the canopy and the axial 
stress in the shroud lines at maximum deceleration. 

5. 1 Hoop Stress in the Canopy 

The maximum pressure on the canopy, calculated from Newtonian 
flow incident at the angle a with respect to the local surface normal: 

p = (pV^) cos^ = 2 S Ivl cos (5.1 

max r max r | (max 

where |^j max i s given by Eq. (4.7). Therefore the hoop stress, in 
accordance with the standard formula for a conical shell of local 
diameter D , is 

p D 3 vl Dcos^a 
max _ r I Imax . 

^h 2 6 sin a 6 sin a ' * i 

where 6 is the effective shell thickness and a = 45° (see Fig. 1). 

We impose the requirement that shall not exceed the allowed stress 

c of the canopy material 


cr, < cr 
h — c 


(5.3) 


Consequently, in view of Eq. (5. 2), 6 must not be less than a minimum 

value 


6 > 6 


min 


= p M 


max 


2 

cos a 


esc 


aj t 


(5.4) 
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The tensile strength or breaking tenacity f of fabrics is usually 

expressed in units of gm/denier , which may be converted to <r 

2 c 
in lbs /ft by the conversion formula 

<r = 2. 95x1 0 4 f p (5. 5) 

c c 


where p is the material density in lbs /ft' 


Density p and f data for the principal plastic under consideration 
here, Nomex, are shown in Fig. 4, obtained from Ref. 4. It is to 
be noted that the tensile strength cr^ at 500°F obtained from this 
figure at f = 3 is cr = 7 . 64x1 0^ lb /ft ^ = 53,000 lb /in ^ , roughly 
one-half that of stainless steel. Other high tenacity plastics, such 
as Dacron, indeed possess tensile strengths comparable with steel 
(Ref. 5). However, because of inferior high temperature properties, 
they will not be further considered here. 


We illustrate the order of magnitude & min to show that the peak 
stress loads can be met with a relatively small effective Nomex thick- 
ness 6 and, consequently, small canopy weight compared with payload 
weight. To express the dependence of 6 on entry velocity 

and ballistic coefficient (3 , we rewrite Eq. (5.4) replacing |v 

from Eq. (4.7) and replacing D by use of the approximation 
fCf. Ea. (3. 2)1 

«- x % ' J 



D = 


>TT §E p ~ Y 


8 W 


^ g E P 


(5.6) 


The result is 


6 . 

min 


< V e sin Ye> 
2 e H M V 


V 


8 W (3 

_e_ 


* § 


E 


2 

cos a esc a 


(5.7) 
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Calculation from Eq. (5.7) with the Nomex datum cr^ = 53, 000 lb/in and 
a = 45° indicates that even for the extreme values of V = l6kft/sec, p = 10 


6 . s 1.4x10 ft 

mm 


(5.8) 


The associated canopy weight is 


c , mm 


u 6 . p 

mm [_c_ 

8 cos a 


~ . 0. 027 lb 


(5.9) 


which is indeed small compared with . A much greater canopy 

weight would be required for the large diameter canopy with low (3 
since W scaleswith (3 according to W ~6 D^~ \/(3 p ^ ~ p ^ ^ 


_ O 

Thus for p = 10 slugs / ft ^ , the minimum canopy weight required for 
meeting the peak hoop load is 0. 27 lb. 


c mm 


5. 2 Axial Stress in Shroud Lines 


The maximum axial stress c r in the shroud lines is associated with 

a 

the maximum deceleration of the combined payload and shroud line 
masses (W^ + W s )/gg . For the assumed configuration in Fig. 1, 
the stress equation is 


W + W 

— 2 5.. 9-^. = A cr sin 60° = it N cr sin 60° 

g E dt max 3 a s s a 


(5.10) 


where the total cross sectional area A is the sum due to N shroud 

s s 

lines of radius R 


In the preceding formula the shroud design parameters N g , R g are 


also implicit in , as given by 


W - tt R N Dp 
s s s s 


(5.11) 


NATIONAL tNiolNttmNS SCIENCE COMPANY 
S-424 


max 


Inasmuch as cr must not exceed the permissible stress level c r in 

SL s 

the shroud material, Eq. (5. 10} represents a limitation on the choice 
of the parameters N g , R g . For the decelerations 
obtained in Eq. (4. 7) it is found that <r < <r can be readily achieved 

3 . S 

with shroud weight W g <<W^ and, therefore, to a high degree of 

N is governed by the constraint 


|dV / dt | 


approximation the choice of R g , 


-> vv 

N R 2 > — 2_ 
s s ~ it g E 


dV 

dt 


max 


o 

esc 60 
cr 

s 


(5.12) 


To estimate the order of magnitude of N g R g , we rewrite Eq. (5. 12) 
in the form 


N 

s 


R 



S ~ IT 



esc 60 


2.95x10 fp 


(5.13) 


and substitute Nomex data at 
result is 


V 

e 


1 6 kft/ sec as in Section 5.1. The 


(N R 2 ) ss 1. 9 x 10" 5 ft 2 
s s 

mm 


(5.14) 


The associated minimal shroud line weights in Eq. 


_ AO VU 

— . 1U 


s , mm 

D = 89 ft . 
max 


tor 


mm 


_ L 1 tf. 

U • J. C O.J 


w 


s , mm 


(5. 1 1) range between 

— A 1 f~\ -v« 

— • TIV J.W A^JL 


It is apparent from the results of this section that the weight restriction 

W + W < W imposes no difficulty in designing to meet maximum 
s c ~ p a o 

dynamic stresses with the M fine-dimension :i criteria of Eqs. (5.8, 5.14). 
Indeed the weight limitation will not be exceeded even with linear dimen- 
sions R , 6 , ten times greater than those obtained here. Greater 
dimensions not only increase the structural safety, but simplify fabrica- 
tion problems not considered within the scope of the present investigation 
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However, as we shall see in the next section, heat transfer and 

material temperature limits impose an upper limit magnitude on 

R and 6 and on the design feasibility consistent with the weight 

restriction W + W < W 
s c ~ p 
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6. SIMPLIFIED HEAT TRANSFER 


The most severe heat load conditions during entry occur, in general, 
along the leading edge of the canopy and the windward element of the 
shroud lines. The heat flux in these critical locations may be due 
either to free molecule flow (fm) regime at extreme altitude, transi- 
tional regime (tr) at intermediate altitudes, or stagnation flow heat 
transfer associated with boundary layer flow ( b 1 ) in the continuum 
regime of the lower altitudes. 

6.1 Transitional Flow 


The various regimes of interest are determined here by assuming that 
transition obtains when the density -dependent mean molecular free path 
X(p) becomes comparable with the dimension R of leading edge or 
shroud radii. In particular we assume that the density range of 
the transitional regime extends over a decade of mean free paths as 
defined by 


— <^M <V To (6.1) 

vTo R 

with mid-transition density p = p obtained at X = R . 

In view of the inverse dependence of X upon the local density p(z) we 
determine from Eq. (6. 1) the relation between mid -transition density 
P - P t and the system dimension of interest R , utilizing standard 
air data at the earth's surface (Ref. 1), 

- 3 3 - 7 

P q = 2.38x10 slugs/ft X Q = 2.18x10 ft (6.2) 


which yield 
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(6.3) 


p* = P X .. R' 1 = 5 . 1 9x1 o" 1 ° R' 1 
r tr r o o 


Accordingly, the densities which define the limits of free molecule 
and continuum flows, in view of Eq. (6.1) are related to R by 


free - molecule 


continuum 


p < p = 1. 63x10" 10 R _1 


p > p = 1 . 63x1 0 R' 1 


(6.4) 


Logarithmic plots of the limiting densities p , p£ m > versus R 

are shown in Fig. 5. The reciprocal density is plotted as ordinate 

in Fig. 5 so that for exponential atmospheres of the type given by Eq. (4. 1) 

the altitude increases in proportion to the ordinate, as indicated in 

the figure. Numerical data for the VM-3 and VM-7 atmospheres 

implies that over the range of "fine" characteristic dimensions 
- 1 - 4 

10 ft , transition may occur from 400 kft down to 

150 kft as R is made smaller. 

6.2 Free Molecule Heat Transfer 


For densities p < p£ m the convective heat transfer is calculated 
from the expression 


„ . . PV (d) 

Q fm^ = ^Tj 


16. 5 ) 


corresponding to diffuse reflection of the free molecule flow incident 

perpendicularly upon a surface having an accommodation factor of 

unity (Ref. 6. p. 43). For the exponential atmosphere, the maximum 

free molecule heat flux Q, is obtained by differentiation with 

fm,max 

respect to density 


^^fm n 4 -j 

— = 0 = d? lo S ° 


fm 


-1 , 

P - 3 c 


(6.6) 
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which is satisfied at the density 


j 2 p sin Y e 

P ~ - P Q, fm " 3 H m 


or altitude 


z 


3 P*, 

rj , M r M 

Z Q, fm “ M 2 p sin y 


(6.7) 


( 6 . 8 ) 


where substitution has been made for c from Eq. (4.4). 

e 

The variation of with p obtained from Eq. (6.7) is also shown 

in Fig. 5. The applicability of the associated expression for the maxi- 
mum heat transfer 


Q 


fm, max 


p sin y 

e 1 e 

3 e H J 
M 


(6.9) 


is restricted to design dimensions R for which p, (R) > p (P) . 

This restriction is satisfied by the smaller values of p — * 10” 3 even 

for relatively large value of R , as shown in Fig. 5. Furthermore, 

it is seen that even for the largest value of p = 10 * the maximum 

value of R permissible in order to assure free molecule conditions 

-4 

at peak heating is achieved at R 5x10 ft , which is appreciably 
greater than the minimum required to satisfy the stress as given by 
Eqs. (5. 8, 5.14). 

6 . 3 Continuum Heat Transfer 

For densities p > p, , the convective heat transfer is calculated 
r r bl 

from an expression applicable to high-speed stagnation flow (Ref. 1, 
pp. 96-97) normal to a "cold" cylindrical body of sectional radius R 
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( 6 . 10 ) 



In the form given here the temperature-dependent properties have 

been accounted for by the numerical factor 645 consistent with the 

2 

units V ft/sec, R ft and Q, , BTU/ft sec . We note that, in con- 

bi 

trast with the referenced heat flux formula, the stagnation enthalpy 
at the wall has been assumed negligible compared with the high-speed 
free-stream stagnation enthalpy. Furthermore, the numerical factor 
645 appropriate for a cylinder, rather than 865 appropriate for a 
sphere, has been introduced here, consistent with the so-called 
Mangier transformation (Ref. 7) which interrelates the boundary 
layer effects in the two geometric configurations. 

By substitution of density and trajectory relations into Eq. (6. 10) 
and differentiation with respect to p 


dQ hl _1 

— =° = ^ l °gQbi = V- 3 - l5 S 


(6.U) 


we obtain the maximum heat flux 


Q 


( v e ) 3 ’ y 


= (6.3=ep/ 5 


bi-max 7 ' e 1 o’ 


( 6 . 12 ) 


at the density 


1 _ P sin Y e 

P ‘b.3c o = P Q,bl = 3.15 H. . 

e M 


(6.13) 
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which obtains at altitude 


z 


6.3 H p 

tt 1 M r M 

S Q, bl " M g 2 p sin Y 


(6.14) 


We note that, in contrast with the result for Q , in Eq. (6.9), 

fm, max n 

the result for Q n , indicates an explicit dependence on the leading 

bl, max 

edge radius of curvature R . Furthermore, ^ formally 

smaller than p,_ , and therefore z_ , . is formally greater than 
^Q, fm Q, bl 7 6 

z_ . in accordance with 
Q, fm 


Q, bl 


z , + H 

Q, fm M 


log 2. 1 


(6.15) 


The dependence of p^ ^ and z q on P is also indicated in Fig. 5. 

The fact that p^ ,, occurs at higher altitudes than p_ , implies 
Q,bl 6 r Q, im 

that if R is made sufficiently small to assure peak heating before 
transition to continuum, the subsequent continuum heat transfer will 
remain well below the peak value described by Eq. (6. 12). 


6.4 Transitional Heat Transfer 


Regarding the transitional regime, no formulas are available at 
present for the heat flux Q. . However, as noted in Ref. 6, p. 59, 
because of the continuous behavior of heat transfer on decreasing 
Reynolds number in going from the boundary layer flow to the free 
molecule flow, the results for Q in engineering practice may be 
interpolated between those of and C>£ m • We show in the follow- 
ing paragraph that the ratio Q f is of the order of unity through- 

out the transitional regime, and utilize this fact in postulating an 
interpolation formula for the transitional heat flux Q 

tr 
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The ratio Q^/Q^ obtained from Eqs. (6.5) and (6.10) is 


Q 


fm 


(P) - 


10 12 v/p p R 
r o 


*bl 


1290 J V/10 


(v/i o 4 ) 0- 


15 


( 6 . 16 ) 


At the mid-transition density p = p given by Eq. (6.3) we have in 
the approximation that exp . 15 c g p — 1 


m 


io 12 p vr 

o o 


P = P tr 1290 J 


(v/10 4 ) 0 ' 1 


(6.17) 


an expression independent of R . 

A 

Numerical substitutions with a representative entrv speed V = 10** ft/sec 

e 

yield: 


fm 

\l 


~ 1 . 10 


at 


P = P 


tr 


V = 1 0’ ft 

e 


(6.18) 


Within the transitional regime defined in Eq. ( 6 . 4) the heat flux ratio 
for V =10 ft/sec varies in the range 


62 = 


10 


10 

174 


Q 


< 


fm 


Q 


< 1 . 10 ( 10 ) 


1/4 


1.97 


bl 


(6.19) 


Finally since the ratio in Eq. (6. 17) is proportional to V ^ the 

e 

limits indicated in Eq. (6.19) may change as V ranges over the 

3 3 e ° 

interval 6x10 ft/sec < < 16x10 ft/sec . The extreme limits of 

the ratio are obtained by changing 1.96 to 2.12 as 

V -* 6000 ft/sec, and 0.62 to 0. 56 as V —>16,000 ^/ 360 . 
e e 


Based on the limited variation of the ratio Q^/Q^ through the transi- 
tional regime, we employ, where necessary, an interpolation formula 
for calculating Q as described below with the aid of Fig. 6. We 
consider the entry system with given p and leading edge dimension R. 
The given p determines the p^ and ^ lines introduced in 
the figure. These lines intersect the lines representing the transition 
density limits as shown in the figure, with point A at the intersection 
of and Pq ^ , and point B at the intersection of p- b ^ and 

p , . . The radii R at which A and B are obtained is designated 

V-6 y D X 

by R^ , Rg . Now if the given R is greater than R fi , then transition 

precedes the peak heat flux (at density p^ ^ ) an< ^ the magnitude of this 

heat transfer is determined unambiguously as Q, , = Q K1 (P, R) 

given by Eq. (6. 12). Similarly, if R is less than R^ the peak heat 

flux is obtained before transition at the density p^ at the peak 

value Q fm,max = “im.J 1 given by Eq. (6.7). If the given R is 
between R and R we employ the following interpolation to calcu- 

n -LJ 

late the transitional peak heat flux 



( 6 . 20 ) 


This interpolation form is constructed so as to satisfy the dominant 
regimes approached as R — ► R^ , and R Rg > Q tr —* • 

Furthermore, the use of a logarithm rather than linear interpolation in 
R suppresses the relative importance of this parameter, consistent 
with its role in the logarithmic plots of Fig. 5. 
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R 


R 


A 


R 


B 


Figure 6 

Sketch of interpolation zone for postulated transition 
heat flux formula Eq. 6. 20. Interpolation required if, for 
specified R and j3, the characteristic dimension R 


is intermediate between R 


obtained from 


P Q.£m< R A>] and R B [ 


obtained from P b i^) - pQ,bl (R B*] ' 


PfmW = 

Q, bl' 


PA- 


r 
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7. FEASIBILITY LIMITS 


As we have shown in Section 5, if the tensile strength of the chute 
material is of the order of 50,000 to 100, 000 psi, then the structural 
weight requirements (W + W ) for supporting the maximum decelera- 
tion loads are small compared with the payload weight . However, 

the temperature dependence of tensile strength imposes a limit on the 
allowed material temperature, which must therefore remain below a 
specified maximum value during entry. We therefore consider 

the feasibility limits which relate the allowed temperature T to 

entry heating. In particular we determine the relation between T 
and the entry velocity V g which, as shown in Section 6, governs the 
maximum heat loads imposed on the vulnerable elements of the chute 
design; i.e. , on the leading edge of the canopy and on the windward 
elements of the shroud lines. Thus we determine the feasibility 
limits for two essentially different deceleration system designs. 

The first of these, having "fine" characteristic dimension R 
(canopy leading edge radius or shroud line radius) will experience 
maximum heat transfer in the free molecule flow of high-altitude 
regime, so that when continuum flow is established at the lower alti- 
tude, convection heating is well below its possible peak magnitude. 

In this system the relation between T and V is found to be 

independent of the (small) deceleration system weight. The second 
system envisages "blunt" dimensions R which incur transitional 

and/or continuum flow heat transfer at peak magnitude. Because 

1/2 

this peak magnitude depends inversely on R , it can be mitigated 

only by increasing R , and therefore, by increasing the deceleration 

system weight. Thus the relation determined between T and V 

max e 

for the "blunt" design depends significantly on the deceleration system 

weight. It is found then that a dominant system parameter which 

appears in the feasibility limit relation connecting T and V 

° max e 

is the weight of the "torus" forming the "blunted" leading edge of 
the canopy. 
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7 . 1 Feasibility Limits with Maximum Heating Q£ m max 

7.1.1 Stress Constraints 

It is seen from Fig. 5 that even with the largest value |3 = 10 *, the 
peak heating will occur in the free molecule flow regime, provided 

- 4 

the "fine" dimension R does not exceed an upper limit 5x10 ft. 

riiaX 

Lower limits on these dimensions arise from strength requirements as 
expressed by Eqs. (5. 8, 5. 14) for Nomex parachutes under the most 
severe entry conditions under consideration. Therefore the feasible 
design dimension by which peak heat transfer is incurred in the free 
molecule flow regime is constrained in the first place by the restriction 

2x5x1 O' 4 ft > 6 > 6 . =1.4xl0 _5 ft (7.1) 

— — mm 


5x1 0" 4 ft > R g > V - ’ffi— (7.2) 

S 

valid for all S < 10 ^ and V <16 kft/sec. 
r — e — 

As an example of satisfactory design data from the stress viewpoint 

for (3 = 10 ^ we may select 6 = .001 in; R = .0036 in for N = 300. 

The dimension D is, then, D = 76.5 ft. The corresponding canopy 

weight W = 23.3 lb and shroud weight W = . 6 lb total less than 
c s 

the maximum permissible chute weight of 50 lb. These results hold 

“ 3 

for Nomex with c r ss 50,000 psi and density p c = 86.3 lb/ft . 

It is to be noted that the stress conditions for metals of greater tensile 
strength have smaller lower iimits than those for Nomex in Eqs. (7.1) 
and (7.2), so that despite the greater metallic densities, the strength 
requirements can be achieved with + W g « W . Furthermore, 
the high temperature strength of metals enables extension of feasibility 
limits with increasing entry velocity beyond those determined below 
for Nomex. 
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7.1.2 Temperature Constraints 


Next we estimate the entry velocity limitation arising from peak heat 
loads. We rewrite Eq. (6.7) in the form 


um , max 


P V e Sin Y e 
3 eH M J 


(7.3) 


and consider its effect on the canopy leading edge and on the shroud 
lines . 


In applying Q„ to the leading edge we consider that, because of 

rL 1 b fm, max 

the effective bluntness 6/2 of this element, the mean heat load is 

— Q ir D ( it — I . The averaging performed here allows for 

2 fm, max 2 

the inclination of the blunted leading edge surface elements to the 

direction of the incident flow. Gray-body re-radiation from the sur- 

4 6 

faces affected by the incident flow is given bycr'eT tt D (ir -^r) where 

T has been assumed uniform throughout the leading edge cross section 

by virtue of its small dimensions and good thermal conductivity. 

Therefore under the assumption of quasi- equilibrium between mean 

heat influx and gray -body re - radiation, we obtain the following equation 

for the maximum temperature T in the leading edge 

c max ° 


cr'e T 4 = j Q f 

max 2 fm, max 


(7.4) 


Analogous considerations applied to the shroud lines indicate 


cr’e T 4 = ^ Q, cos 2 60 

max 4 fm, max 


(7.5) 


Therefore the higher temperature calculated from these equations, i.e., 
that given by Eq, (7.4) for the leading edge, governs the feasibility 
limit in free molecule flow. 
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Substituting from Eq. (7.4) into Eq. (7.3) we have, upon simplifying, 

the following relation between T and V : 

max e 


T 

max 


(3 sin y e 

6<r' € e H J 
M 



v 3/4 


(7.6) 


Substitution of numerical values 
o-' = . 476x1 0 _1 2 BTU/ft 2 sec (deg R) 4 
sin Y e = • 259 e = 2. 72 


J = 778 ft-lb/BTU 

H. , = 50. 5x 10 3 ft 
M 


results in 


1 / 4 T 

max 


5.37 p 1 /4 V 3/4 
r e 


(7.6') 


Curves of T max versus V e obtained from Eq. (7. 6') for several 

values of (3 are shown in Fig. 7. A tabulation of maximum entry 

velocities versus (3 for Nomex with T = 500°F with assumed 

max 

value of e = . 92 

Table of V for Nomex 

e, max 

(3 slugs /ft 2 1 0 1 

V ft/sec 2090 

e , max 


10" 2 10" 3 
4510 9700 


-3 

The minimum value of the chute weight for the p = 10 design, 
obtainable from Eqs . (5.8, 5.11) is 


W . + W 

c , min s 


min 


= 0. 27 + 0.46 = 0. 73 lb 


(7.7) 
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7.2 Feasibility Limit with Maximum Heating max 

As noted above the required deceleration system weight is now dependent 
on preventing excessive temperature in the blunted leading edge of the 
canopy. The "blunting" of the shroud lines need not incur much addi- 
tional weightj since the shroud lines are favorably inclined to the 
direction of flow [Cf. Eqs. {7.4, 7.5) for the effect of inclination]. 

Accordingly, we introduce the parameter representing the weight of 
the deceleration system, as consisting essentially of the ring weight: 

P = W /W (7.8) 

r p ' ' 


in terms of which we may conveniently express other design dimensions, 
such as diameter D , and the ring radius : 


D = 


V 


8 W p (1+p) 

* Sr P 


(7.9) 



(7.10) 


The latter dimension is of dominant importance in the maximum 
stagnation heat flux Eq. (6. 1 2). 


We now calculate the maximum temperature with qua si -equilibrium 
between continuum heat in flux and thermal re-radiation [Cf. Eq. (7.4)] 


cr 1 e T 


4 

max 



bl, max 



(7.11) 
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where Eqs. (6. 12, 6. 13) have been used and allowance has been made 
for re-radiation along the whole cross sectional circumference of the 
ring. Further, we may express the latter equation in the form 


4 3 15 

o-'eT* = F (fl, p) V 
max r e 


(7. 12) 


where F is a function of the system design parameters p , p . 

, 1/8 

l A i /a V/^/r/ P t 

F ((3 


/ 7 \ 1/8 / 1 

, 161 (S\ l/8 V^( P r ) / Sin \ 

TT^w "'\g E v V 3 - 15 eH M p o/ 


1/2 /l+p\ I/8 3/8 
M Po/ (p 2 ) ' 

(7.13) 


Therefore, the feasibility limit equation, Eq. (7.12), rewritten in the 
form 


1/4 T _ p 1 / 4 y 0. 7875 

e max e 


implies a weak dependence of on p implicit in F 


(7.14) 


1/4 


Limit curves calculated from Eq. (7. 14) for Nomex density = 86.3 lb/ft“ 
are shown in Fig. 8 for several values of p and (3 . Continuum con- 
ditions are satisfied if >VT0 \ (p) which, by use of Eq. (7.10), may 
be converted to an expression relating p and p 


4/5 ? 1/5 

J 10 \> P o eH M\ ( P r \ / 1 +p\ X 

V ) W W p/ 1 p 2 J 


/ 5 


(7.15) 


The special form of Eqs. (7. 14) and (7. 15) with Nomex data under the 
entry conditions of interest is 
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Feasibility limit curves for entry heat pulse in continuum and transition flow 





y 0. 7875 


(7.1 4') 
(7.15') 


From these expressions we obtain the limit curves showing as solid 

1 /4 

curves the variation of « T versus V in Fig. 8. In the same 

max e 

figure, the transitional flow limits are shown in dashed curves. We 
note that subject to the restriction of the weight parameter p < 1 

continuum can be obtained, in view of Eq. (7. 15'), only if (3 > 1 . 9 x 

- 2 2 - 1 
10 slugs /ft . Thus, the results in Fig. 8 show the (3 = 10 family 

-2 

in continuum flow for the indicated values of p < 1 ; the (3 = 10 in 

-2 

transition flow of indicated values of p < 1 . For (3 = 2x10 continuum 
is obtained at p = 1 , but for p = . 5 , transitional flow will obtain at 
maximum entry heating. 

It is to be noted that the (3 = 10 ^ does not appear in Fig. 8, because 
for the Nomex system all possible configurations subject to p < 1 will 
experience peak heat load in free molecule flow. Thus, the maximum p 
of interest in the boundary layer regime is about 2x10 . Accordingly, 

as shown in Fig. 8, the maximum feasible entry velocity in "laminar" 
peak heating is 6 330 ft/sec. 


Finally, we note from Figs. 7 and 8 that if entry at small speeds V 

- 1 e 

is to be achieved with the higher (3 system ((3 — ♦ 10 ) the blunted 

design, rather than the fm design, becomes feasible, but for high- 

-3 

speed entry the fm design with (3 — * 1 0 governs the feasibility 
limit. 
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8. DISCUSSION OF RESULTS 


The feasibility limit of high-speed entry with all-Nomex parachutes 
is obtained from the curves in Figs. 7 and 8 by locating the value of 

V at the ordinate 
e 

e 1/4 T =;(.92) 1/4 (460 + 500) - 940°R 

max — 

The value of € ^ 0. 92 used here is presumably obtainable ' by appro- 
priate coating or processing of Nomex. It is thus seen that to achieve 

high entry velocities with Nomex, it is necessary to have the lowest 

-3 

possible ballistic coefficient (3—10 , with "sharp" leading edges 

(windward elements) which experience peak heat load in the free 
molecule regime. The feasible entry velocity from Fig. 7 is then 

V css 9 , 700 ft/sec , well above the feasible maximum value obtained 
e 

from Fig. 8 in laminar (continuum) heat transfer V — 6, 730 ft/sec , 

- 2 e 

(3 = 2x10 . In the latter (bl) design, the parachute weight is 50 lbs 

(p = 1) , incurred mostly in the weight of the toroidal "blunt" leading 

edge of the canopy. On the other hand, the former (fm) design 

requires a minimum weight of only [Cf. Eq. (7.7)] 

W . +W . as (0. 27 + 0.46) = 0. 73 lb (8.1) 

c , rmn s , mm 

or (i = 0^5" • The weight and maximum velocity advantage of the fm 
versus the bl parachute design can be ascertained by comparing the 
temperature maxima obtainable under these flow conditions, using 
Eqs. (7.6', 7.14'): 


* "Study of Expandable, Terminal Decelerators for Mars Atmosphere 
Entry, " Goodyear Aerospace Corporation, Akron, Ohio, October 3, 
1966 
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( 8 . 2 ) 


max, fm 
T 

max, bl 



y -. 0375 


This ratio decreases as the chute weight associated with the mass 
parameter p decreases, as (3 decreases, and as the entry velocity 
increases . 


The relatively small minimal weights for the fm design with Nomex, 

W . + W « 50 suggest the possibility of using denser, high 

c,mm s,min && } 6 s 

temperature high-strength metals for further extending the feasible 
entry velocity within the parachute weight restriction of 50 lb. Con- 
sider, for example, stainless steel (SS 301) which retains a tensile 
strength of about 60, 000 psi at the high temperature of 1000°F. For 
this temperature we obtain from Fig. 7 the following table of feasible 
entry velocities 


P 




10 


-3 


V 

e 


3660 7890 17,000 


To assure meeting both stress requirements and fm conditions at 
maximum heat transfer, it is sufficient to keep the same design 
dimensions as calculated above for Nomex. Accordingly, for example, 

-3 

for (3 = 10 the minimal weight requirement obtained by scaling the 
result in Eq. (8. 2) with the steel-to-Nomex density ratio of 5. 7 is 


W . + W . =1.5 + 2. 7 

c , min s , mm 


(8.3) 


Implicit in the minimal weight W , is a minimal canopy thickness 

^ C ; I XT X XI 

^min ~ It , in view of the scaling formula Eq. (5.7). From 

the fabrication viewpoint such a minimal thickness may indeed be 
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problematic. We note, however, that the canopy weight may be 

increased from its minimal value of 1 . 5 lb in Eq. (8. 3) to 47.3 

without exceeding the parachute weight restriction of 50 lb. With 

the indicated increase, the effective thickness 6 will also increase 

-5 -4 

to a practicable value of 6 a (47 .3/1.5) 6 cs 4. 4x1 0 ft = 5 . 3x1 0 in. 

Finally, we note that the T as discussed in Section 6 generallv 

occurs at a higher altitude (i.e., earlier) than the peak stress. Cool- 
ing in the interval between peak heat load and peak stress permits an 
increase in structural strength above that employed in the calculations. 
From this viewpoint the calculations are somewhat conservative. How- 
ever, since T depends on the fourth root of the heat flux, inclusion 

max 

of the strength recovery effect in the calculation does not appear 
warranted. 
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9. REMARKS ON SHOCK IMPINGEMENT 

We consider the shock wave which may impinge on the parachute due 
to the bow shock formed by the primary body. 


The primary body diameter 
with D (see Appendix) 


D is assumed to be small compared 


D =* V 10‘ 3 D = . 0316 D (9.1) 

P 


The profile of the (hypersonic) shock wave, formed at altitudes where 
» X. , will be assumed in the approximate parabolic form, given 
by blast wave theory (Ref. 8) 



(9.2) 


where 


x = axial distance measured from nose 

y = lateral (radial) distance from axis 

2 1 = spherical nose drag coefficient 
P 

In the "exit" plane of the conical canopy the radius of the shock wave is 


0. 9 (D D) 1/2 = .16D 

P 2 


D 

2V2 


(9.3) 


Therefore the axisymmetric shock front will lie well within exit hole 
of the canopy, without giving rise to an impingement problem. 


41 


K i (H. f « •' — N K < A 1 

» IUINML 


ENGINEERING SCiENCc. uOmfanV 


S-424 


A portion of the shroud lines near their attachment to the primary- 

body will in general be submerged under the bow shock wave 

developing at the lower altitudes, and in the "entropy layer" close 

to the body axis. From the geometry of the shock profile Eq. (9. 2) 

impingement will occur as a "weak shock" at distance of radial 

distance of about D from the axis, and an axial distance of 3 D 

P P 

downstream of the nose. Although the temperature is high in the 

wake of the shock and entropy layer, the density is in fact smaller 
than ambient density. Hence the fm design with its fine shroud 
lines may remain in free molecule flow well below the altitudes 
where the primary shock has developed, and its heating rate is not 
enhanced under these conditions. Analogous considerations indicate 
that the "blunt" design with fewer shroud lines will not experience 
more severe heating in the "weak shock" and entropy layer than 
due to that of the ambient flow. Finally to assure that the fine 
shroud lines of the fm design will not be submerged under the 
shock wave profile given by Eq. (9. 2), they may be terminated on 
a few short "blunt" lines piercing the bow shock profiles and con- 
tinuing the connection to the primary body. Details of such a com- 
posite attachment design would require further study. 
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10. CONCLUSIONS 


The following conclusions are made on the basis of the present investi- 
gation. 


la. The axisymmetric Nomex parachute will survive laminar 


entry heating pulse up to an entry velocity V 


= 6330 
max 7 

ft/sec with a ballistic coefficient of (3 = 2x10^ slugs/ft“ 
corresponding to maximum deceleration system weight of 
50 lb. Most of this system weight is due to the ’’blunt” 
toroidal ring at the leading edge of the canopy. 


lb. With the weight restriction as above, an entry system 

designed to undergo entry heat pulse under free molecule 
flow conditions can survive entry heating up to entry 
velocities V_ as listed below for two feasible high- 
temperature materials. 


Maximum V g with Heating Pulse in Free Molecule Flow 


2 

(3 slugs /ft 

V ft/sec (Nomex at 500°F) 
e 

V e ft/sec (Steel at 1000°F) 


lo " 1 

10" 2 

io “ 3 

2090 

4510 

9700 

3660 

7890 

17,000 


2. With the weight and entry velocity restrictions as above, 
the heated parachute can survive the peak dynamic loads 
arising in all flow regimes considered, 

3. The Mars atmosphere entry conditions are favorable to 
the possible utilization of "sharp" leading edge decelera- 
tion systems which enable critical entry heating in the 
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free molecule flow regime. As a consequence, all-metal 
deceleration systems, possessing strength at high tempera- 
ture, provide an approach to the extension of entry velocity 
limits, without incurring excessive weight penalties. 


44 


SCIENCE COM PAN V 


S-424 


11. ACKNOWLEDGMENTS 


The writer is pleased to acknowledge helpful discussions with 
E. Laumann and B. Dayman of JPL and with G. Bergman and 
J. Grzesik of NESCO. 


45 

NATIONAL ENGINEERING SCiENCE COMPANY 


S-4Z4 


12. REFERENCES 


1. Martin, J. J. , Atmospheric Reentry, Prentice-Hall, Inc., 
Englewood Cliffs, N. J. , (1966). 

2. Loh, W. H. T. , "Ballistic Reentry at Small Angles of Inclination, " 
J. American Rocket Society, Vol. 32, (1962), p. 718. 

3. Gazley, C. , Jr. , "Heat Transfer Aspects of the Atmospheric 
Reentry of Long Range Ballistic Missiles, " Report R-273 
(Santa Monica, California: The RAND Corporation, August 1, 
1954). 

4. E. I. Du Pont de Nemours Co. - Bulletin N-201, "Properties of 
Nomex, " October 1966. 

5. E. I. Du Pont de Nemours Co. - Bulletin D-173, "The Tenacity 
and Elongation of Dacron, " January 1964. 

6. Probstein, R. F. , "Heat Transfer in Rarefied Gas Flow, " 

Theory and Fundamental Research in Heat Transfer, The 
MacMillan Company, (1963), Ed. J. A. Clark. 

7. Schlichting, H. , Boundary Layer Theory , McGraw-Hill Book Co. , 
Inc., New York, (19b0), p. 190. 

8. Lykoudis, P. S. , "Ionization Trails, " Proceedings of the 1961 
Heat Transfer and Fluid Mechanics Institute, Stanford University 
Press, Stanford, California, (1961). 


46 


S- 424 


PRECEDING PAGE BLANK NOT FILMED. 


APPENDIX 




S-424 


preceding page blank not- filmed. 

PRECEDING PAGE BLANK NOT FILMED. 

Exhibit No. 1 

' RFQ BE- 38465 

21 March 1967 


A STUDY OF THE USE OF A 
PARACHUTE ENTRY SYSTEM FOR MARS LAUDER 


Introduction 


The relatively tenuous atmospheric conditions expected to exist at Mars 
present a significant engineering challenge in determining the optimum design 
for a s mall survival scientific entry probe. Studies have been conducted of 
blunt body aerodynamic decelerators, single and multistage subsonic and super- 
sonic parachutes, guided retro rocket landing systems and lifting entry. Each 
of these systems have certain drawbacks. 

An alternate entry system has been proposed; i.e., a payload suspended 
from a large diameter ring parachute which was opened outside the sensible 
atmosphere of Mars. This system could conceivably have a lower impact 
velocity than the blunt body devices, not require the Mach and/or altitude 
sensors required for parachutes and would avoid the opening shock problems 
and would not require the altitude and directional control systems required 
for retro and lifting entry. On the other hand several significant questions 
can be raised concerning the feasibility of such a system. Two questions in 
particular require answers before any detailed feasibility study can be con- 
sidered. 

The two questions requiring answers are: 

1. Assuming state-of-the-art materials are used, can a parachute sur- 
vive the entry heating pulse? 

2. Can this heated parachute survive the peak dynamic pressure pulse 
while carrying a meaningful payload? 

A first order study shall be made to answer these two questions. 

Assumptions 

Certain assumptions must be made in order to limit the scope of the study. 
A number of basic assumptions* are listed below. 

1. Parachute configurations 

a. Conical, flat annular with central open area equal to l/2 overall 

frontal area (A „) 
o 

b. Shroud line length equal to chute maximum diameter 


*Some of these assumptions are arbitrary. If it is shown that any of them are 
unnecessarily restrictive or avoidable, they may be changed upon request. 
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c. C D = l/2 (based upon overall frontal area) 

d. Weight may not exceed weight of payload 


2 . Payload configuration 


a, 

b. 


Heir.i spliss? 


-cylinder with frontal area (A ) = .001 A„ 
* p 0 


Payload weight of 50 lbs. 


3 . Loads 

a. Equal distribution over canopy 

b. Uniform shroud line loading 

c. Parachute opening shock « q TTigy loading 

4. Heating 

a. Critical peak heating is laminar, convective to cylindrical chute 
leading edge or to shroud line at chute attachment point 

b. Heating on rest of chute is not critical but assume that chute 
panel temperatures of 250° F exist at q m&x 

c. Radiation is outward from grey body to black space 

5 . Materials 

a. Payload container made of non-ablative material 

b. Chute and shroud lines of Dacron and Nomex 

c. Do not be concerned about minimum gages of materials 

Study Parameters 


The effects of variation of several mission parameters should be investi- 
gated in the study. These are: 

1. Ballistic coefficients ( slugs /ft ) between 0.001 and 0.10 

2. Trajectory angle at 800,000 ft is -15° only 

3. Velocity at 800,000 ft between 6,000 and 16,000 ft/sec 

4. Mars atmosphere described by VM-3 and VM-7 engineering models 
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Technical Resources Available 


The JPL will support this study by making available trajectory calculations 
which describe the dynamic pressure and the ideal gas stagnation point heating 
rate histories. Various other information can be made available such as prior 
parachute and drogue studies and characteristics of fabrics at elevated 
temperatures. 


Study Results 

1. Describe the assumptions that were required. 

2. Determine under what parametric combinations that the mission is 
feasible. 

3. For feasible missions describe (giving estimates of the certainty 
tolerance) : 

a. Parachute weight and material 

b. Leading edge diameters 

c. Shroud line sizes 

d. Critical temperature histories 

e. Other unique characteristics 

4. Recommend alternate considerations which could lead to further study. 


Presentation of Results 


Intermediate oral presentations shall be made at JPL after two weeks and 
four weeks of effort to describe progress and plans to date. At the completion 
of the study, both an oral presentation and a memorandum final report will 
be due at JPL. 
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